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ABSTRACT 

The Cosmic Lens All-Sky Survey (CLASS) is an international collaborative program 
which has obtained high-resolution radio images of over 10000 flat-spectrum radio 
sources in order to create the largest and best studied statistical sample of radio- 
loud gravitationally lensed systems. With this survey, combined with detailed studies 
of the lenses found therein, constraints can be placed on the expansion rate, mat¬ 
ter density, and dark energy (e.g. cosmological constant, quintessence) content of the 
Universe that are complementary to and independent of those obtained through other 
methods. CLASS is aimed at identifying lenses where multiple images are formed 
from compact flat-spectrum radio sources, which should be easily identifiable in the 
radio maps. Because CLASS is radio-based, dust obscuration in lensing galaxies is 
not a factor, and the relative insensitivity of the instrument to environmental condi¬ 
tions (e.g. weather, “seeing”) leads to nearly uniform sensitivity and resolution over 
the entire survey. In four observing “seasons” from 1994-1999, CLASS has observed 
13783 radio sources with the VLA at 8.4 GHz in its largest “A” configuration (072 
resolution). When combined with the JVAS survey, the CLASS sample contains over 
16,000 images. A complete sample of 11685 sources was observed, selected to have 
a flux density of at least 30 mJy in the GB6 catalogue at 4.85 GHz (spanning the 
declination range 0° ^ (5 < 75° and |6| ^ 10°, excluding the galactic plane) and a 
spectral index a ^ —0.5 between the NVSS at 1.4 GHz and the GB6. A typical 30- 
second CLASS snapshot reached an rms noise level of 0.4 mJy. So far, CLASS has 
found 16 new gravitational lens systems, and the JVAS/CLASS survey contains a total 
of 22 lenses. The follow-up of a small number of candidates using the VLA, MER¬ 
LIN, the VLBA, and optical telescopes is still underway. In this paper, we present 
a summary of the CLASS observations, the JVAS/CLASS sample, and statistics on 
sub-samples of the survey. A companion paper presents the lens candidate selection 
and in a third paper the implications for cosmology are discussed. The source cat- 
alogues from the JVAS/CLASS project described in this paper are available from 
attp://www.jb.man.ac.uk/research/gravlens/ . 
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1 INTRODUCTION 

The Cosmic Lens All-Sky Survey (CLASS) is an interna¬ 
tional program to map more than 10,000 radio sources in 
order to create the largest and best studied statistical sam¬ 
ple of radio-loud gravitationally lensed systems. CLASS is 


aimed at identifying lenses where multiple images are formed 
from compact flat-spectrum radio sources. The resulting 
lens configurations should be easily identifiable in the ra¬ 
dio maps. CLASS is most efficient at finding galaxy-mass 
lenses with arcsecond image separations, as these mass dis- 
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tributions dominate the leasing optical depth. In principle, 
CLASS can detect more extreme leased systems with larger 
separations due to clusters of galaxies (Phillips et al. 2001), 
as well as smaller-separation (160-300mas) lenses (Augusto, 
Wilkinson & Browne 1998; Augusto & Wilkinson 2001). Set¬ 
ting limits on gravitational leasing on these smaller and 
larger angular scales is a secondary goal of the survey. 

The Very Large ArrayQ (VLA) was used as the primary 
instrument for the CLASS survey. In its largest “A” con¬ 
figuration, the VLA provides high-quality snapshot images 
with angular resolution of 0.2-0.3 arcseconds at an average 
observing frequency of 8.46 GHz. In the four primary phases 
of CLASS carried out from 1994 - 1999, a total of 13832 tar¬ 
get radio sources were mapped. From among these, hundreds 
of multiple-component sources were identified and followed 
up with higher resolution using MERLIN (50 mas resolu¬ 
tion) and the VLBA (2 mas resolution), and with multi¬ 
wavelength observations using a number of instruments. So 
far, CLASS has discovered 16 new radio-loud gravitational 
lens systems, which when added to the 6 lenses from JVAS 
makes a total of 22 confirmed lens systems in JVAS/CLASS 
(Table |^. Some follow up observations are still underway 
but it is doubtful that more multiple image systems will 
be found in the sample. For example, we have yet to con¬ 
clusively identify the CLASS-2 source CLASS B0827-I-525, 
a 2'.'82 double, as either a ‘dark lens’ or binary radio-loud 
quasar (Koopmans et al. 2000b). The details of this follow¬ 
up program and the verihed lenses will be given in Paper 2. 

There have been previous VLA-based lens surveys: the 
MIT-Green Bank (MG) survey (Burke 1989), and the Jo- 
drell Bank VLA astrometric source survey (JVAS) (Patnaik 
et al. 1992a; Browne et al. 1998; Wilkinson et al. 1998). 
JVAS sources, like CLASS sources, were required to have flat 
radio spectra, but the JVAS survey was restricted to 5-GHz 
radio flux densities ^ 200 mjy, compared to the 25-30 mjy 
limit of CLASS. JVAS is thus effectively a bright sub-sample 
of CLASS. The main difference between JVAS/CLASS and 
the MG survey is that the latter has no spectral index se¬ 
lection, and so also picks up sources dominated by steep- 
spectrum lobe emission. Of the ~ 5000 target sources ob¬ 
served in the two previous surveys, hve multiply imaged lens 
systems were found in the MG survey and six in JVAS, of 
which one was a rediscovery of one of the MG lenses. CLASS 
has more than doubled the number of radio-discovered lens 
systems known, and JVAS/CLASS as a whole constitutes 
the most comprehensive single sample of lenses for use in 
astrophysical applications. In addition to these Northern ce¬ 
lestial hemisphere surveys, the VLA is also being used to 
carry out a gravitational lens survey in the Southern sky 
(Winn et al. 2000). 

One of the most important uses of galaxy-mass lenses is 
in the determination of cosmological parameters. In partic¬ 
ular, time delays measured between the components in mul¬ 
tiple image systems with reliable mass models can directly 
determine the angular-diameter distances to the lens and 
lensed objects, and therefore the Hubble constant Hq. Be¬ 
cause the leasing equations depend on a ratio of distances, 


^ The National Radio Astronomy Observatory is operated by As¬ 
sociated Universities, Inc., under cooperative agreement with the 
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the delays are not as sensitive to the matter density 
and cosmological constant Aq. Of course, the success of the 
method depends upon the presence of measurable variabil¬ 
ity in the lensed source, and the ability to construct a well- 
constrained mass model for the lens, which in turn should 
ideally consist of a single deflector well-fit by a spherical 
or ellipsoidal potential. Compact flat-spectrum lensed radio 
sources are more likely to be variable and are easily recog¬ 
nized as multiply-imaged. Steep spectrum, extended lensed 
sources can provide more constraints on the lensing poten¬ 
tial, but are unlikely to be variable and are hard to distin¬ 
guish from complex unlensed radio lobes. Thus, CLASS is 
aimed at filtering out the latter in preference for the for¬ 
mer. Unfortunately, only a few lenses found in CLASS will 
satisfy these criteria and be suitable as cosmological stan¬ 
dards. The CLASS group has successfully monitored the 
quadruple lens CLASS B1608-I-656 (Fassnacht et al. 1999b; 
Koopmans & Fassnacht 1999, Fassnacht et al. 2002), and 
all three independent time delays between the four compo¬ 
nents have been measured. For the current standard cos¬ 
mology (Hm = 0.3, Ao = 0.7) a Hubble constant Hq = 
61-65 (galaxy centroid positions) ±2 {2a statistical) ±15 
(systematic) kms“^Mpc“^ has been derived (Fassnacht et 
al. 2002). In addition, the double image systems CLASS 
B1600±434 (Koopmans et al. 2000a, Koopmans & de Bruyn 
2000) and JVAS B0218±357 (Biggs et al. 1999) were mon¬ 
itored and time delays measured, giving Hubble constant 
values of 57/l^^ and 691 kms '^Mpc ^ respectively (er¬ 
rors are 2cr statistical errors assuming singular isothermal 
ellipsoid mass models). The mass profile modeling for these 
systems is being further refined, and new results will be pre¬ 
sented in upcoming papers. 

The observed lensing rate, and to a lesser extent the 
distribution of image separations and redshifts of the lenses, 
provides constraints upon the cosmology, particularly the 
differential volume of the universe versus redshift which is 
largely controlled by the matter density and the dark 
energy density Aq and its equation of state parameter w 
ranging between w = —1 for a cosmological constant and 
ui = 0 for non-relativistic matter. A well-defined sample se¬ 
lection is necessary to understand the statistics of the lens¬ 
ing, which is in turn necessary to constrain Aq (Turner, Os- 
triker & Gott 1984; Kochanek 1996). For this purpose it is 
vital that all the lenses in the observed sample must be iden¬ 
tified. Because CLASS is radio-based, complications arising 
from dust obscuration in lensing galaxies do not occur; and 
because it targets fiat-spectrum radio sources with intrinsi¬ 
cally simple radio structures, real lenses are easy to identify. 
Furthermore, the external conditions do not strongly affect 
the observations and the survey was carried out with nearly 
uniform sensitivity and resolution. All these factors make 
it the best sample overall for statistical studies of lenses. 
However, interpretation of the statistical results from any 
observational survey in the context of theoretical models re¬ 
quires the control of selection effects in the data sample. 
In the case of CLASS, this means studying the effects of 
the lensing on the parent sample of the survey, the result¬ 
ing image geometries, the ability of the automatic mapping 
procedure to identify the lensed images, and possible loss of 
lenses from the survey. In addition, the CLASS sources suf¬ 
fer from the drawbacks generic to radio-selected samples of 
having poorly constrained redshift distributions and source 
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Lens 

Session 

Max. Sep. 

# images 

B0128+437 

CLASS-1 

O'.'54 

4 

B0218-I-357 

JVAS 

0'.'33 

2 -1- ring 

MG0414+054 

MG, JVAS 

2" 09 

4 

B0445-I-123 

CLASS-2 

T.'33 

2 

B0631-1-519 

CLASS-1 

T.'ie 

2 

B0712+472 

CLASS-1 

T'27 

4 

B0739+366 

CLASS-2 

0'.'54 

2 

B0850-I-054 

CLASS-3 

O'.'68 

2 

B1030+074 

JVAS 

T.'56 

2 

B1127+385 

CLASS-2 

01'70 

2 

B1152+199 

CLASS-3 

T.'56 

2 

B1359-I-154 

CLASS-3 

T.'65 

6 

B1422-I-231 

JVAS 

T.'28 

4 

B1555+375 

CLASS-2 

O'.'43 

4 

B1600+434 

CLASS-1 

T.'39 

2 

B1608-I-656 

CLASS-1 

2'.'08 

4 

B1933-I-503 

CLASS-1 

T.'17 

4 + 4 -1- 2 

B1938+666 

JVAS 

O'.'93 

4 + ring 

B2045-f265 

CLASS-2 

T.'86 

4 

B2108+213 

CLASS-4 

4'.'60 

2 

B2114-I-022 

JVAS 

2'.'57 

4 

B2319-I-051 

CLASS-2 

T.'36 

2 
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Table 1. JVAS/CLASS Lenses The separation for lenses with more than 2 components is the maximum separation. The number of 
images listed is the number of radio components believed to be lensed images, not counting sub-components seen on VLBI scales. The 
principal references are given for each lens. The thirteen lenses from the CLASS Complete Sample (section 2.2) are shown in boldface. 


luminosity functions, which must be dealt with both statis¬ 
tically and through follow-up optical observations. 

Data from CLASS also provide important astrophysical 
constraints upon the mass distribution within the lensing 
galaxies themselves (e.g. Kochanek 1995). Koopmans & de 
Bruyn (2000) observed the signatures of microlensing in the 
radio light curves of the two images in the system CLASS 
B1600+434, and deduce that a significant fraction of the 
mass in the inner dark-matter halo of the lensing galaxy 
is in the form of compact objects. Rusin & Ma (2001) use 
the absence of faint ’’odd” images in deep radio maps of 
CLASS lenses to place a lower limit on the inner mass pro¬ 
files of lensing galaxies. Rusin & Tegmark (2001) consider 
the frequency of quadruple image systems in CLASS and at¬ 
tempt to derive a consistent set of astrophysical mechanisms 
that explain this statistic as well as other observational con¬ 
straints. Keeton & Madau (2001) and Phillips et al. (2001) 
constrain the concentration of cluster-mass dark matter ha¬ 
los based on the paucity of wide-separation lens systems 
in the CLASS sample. Cohn et al. (2001) use the lensing 
model for CLASS B1933-I-503 to constrain density profiles. 
Rusin et al. (2002a) use high-resolution radio and optical ob¬ 
servations of CLASS B1152-I-199 to model the galaxy mass 
profile. Finally, observations of the six-image system CLASS 
B1359-1-154 was used to place constraints on the mass dis¬ 
tribution within the small group of galaxies responsible for 
the lensing (Rusin et al. 2001b). Astrophysical applications 
of lensing such as these studies are as important as the cos¬ 
mological uses outlined above, and may ultimately provide 
critical clues to the nature and properties of dark matter in 
galaxies and groups. 

The success of the CLASS survey, and the unprece¬ 
dented observing and analysis speed, was made possible 
through automation of the observation scheduling and map¬ 
ping analysis pipeline. Specialized auto-mapping software 
(see § ^) was developed to carry out these tasks. This has 


allowed us to observe an average of one target source per 
minute with the VLA, and then later to map the calibrated 
data at a similar (now much faster) rate. 

In addition to the lens candidates and lensing-related 
applications, the CLASS survey provides a solid body of ra¬ 
dio information for all the objects. Astrometric (35 mas) 
positions for the sample are determined, which extends 
the MERLIN/VLBA phase-reference source network in the 
northern sky to a spacing of one to a few degrees (except 
in the galactic plane). CLASS has also proven to be useful 
in studies of active galactic nuclei. For example, JVAS and 
CLASS have been used to identify samples of high-redshift 
quasars (Snellen et al. 2001, 2002). CLASS observations have 
supplemented lower frequency surveys in identifying young 
radio sources (Snellen et al. 2000). Marcha et al. (2001) used 
CLASS to identify a large sample of low-luminosity blazars. 

The JVAS/CLASS survey is now complete. This is the 
hrst in a series of papers which describe the survey and its re¬ 
sults, and covers the initial observations, data reduction and 
properties of the survey. A second companion paper (Browne 
et al. 2002) will describe the selection and followup of lens 
candidates. A third paper (Chae et al. 2002) presents the 
cosmological results, in particular the derived constraints 
on the matter density and the dark energy density and its 
equation of state. 


2 THE SURVEY SAMPLE 

To achieve our goal of obtaining ~ 20 new simple lenses 
suitable for cosmological measurements, we must start with 
a parent sample of more than 10000 flat-spectrum sources. 
Moreover, to be easily followed up with higher-resolution 
radio interferometers such as Multi-Element Radio-Linked 
Interferometer Network (MERLIN), the NRAO Very Long 
Baseline Array (VLBA) and the European VLBI network 
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(EVN), the sources should have total compact emission of 
at least 20 mjy. Note that in a lensed system this flux density 
can be split into four or more images. 

To obtain compact multiply-imaged components, we 
target sources with relatively flat radio spectra. We use the 
convention S' oc where u is the frequency of observation. 
We classify as “flat-spectrum” those sources with a spectral 
index of a ^ —0.5 between 4.85 GHz and a lower frequency. 
Technically, we start with our 4.85 GHz sample and reject 
“steep-spectrum” sources that have a < —0.5 versus a lower 
frequency, as in some stages of the survey targets were kept 
in the sample when low frequency measurements were un¬ 
available. 

The CLASS observations were taken over five years: 
CLASS-1 in 1994, CLASS-2 in 1995, CLASS-3 in 1998 and 
CLASS-4 in 1999. Table ^ shows the progression of the sur¬ 
vey from JVAS through the end of CLASS. The selection 
for JVAS is described by Patnaik et al. (1992a). The CLASS 
survey selection is complicated, and we describe it here in 
some detail as it affects the use to which the data products 
may be put. We also define and describe a simpler subset 
(the CLASS complete sample). 

Our starting point for source selection in CLASS-1 and 
CLASS-2 was the 87GB version of the Green Bank Survey at 
4.85 GHz (Gregory & Condon 1991). However, the updated 
survey GB6 (Gregory et al. 1996) became available and was 
adopted as the parent sample; in this survey the source posi¬ 
tions and flux densities are better determined than in 87GB. 
Furthermore, in the later stages of the CLASS survey, we 
adopted the NRAO-VLA Sky Survey (NVSS, Condon et al. 
1998) as the low frequency survey to define the spectral se¬ 
lection. We now describe the observation sessions and the 
respective selection criteria in detail. 

2.1 CLASS-1 and CLASS-2 

In the first phase of CLASS, observed with the VLA in 
February-May 1994, 3220 target fields were mapped. The 
sources were selected from the 87GB catalogue (Gregory & 
Condon 1991) which contains 54579 sources with 4.85 GHz 
flux densities of 25 mJy or greater in the declination range 
0° ^ (5 ^ 75°. We imposed a galactic latitude limit of 
\b\ > 10°. Spectral selection was based on two-point spectral 
indices versus the 325 MHz Westerbork Northern Sky Sur¬ 
vey (WENSS, Rengelink et al. 1997), the 365 MHz Texas 
Survey (Douglas et al. 1980), and the 1.4 GHz Green Bank 
survey (White & Becker 1992). When available we used the 
WENSS since the flux densities are reliable and sources 
have a positional accuracy of ~ 5 arcsecond. At the time 
of CLASS-1, only a few “mini-survey” areas were complete, 
but by the time of CLASS-2, WENSS was complete above 
S > 30°. In CLASS-1 and CLASS-2, sources outside the 
area covered by WENSS at that time were selected using 
the Texas 365 MHz catalogue. 

WENSS has a flux density limit of 15 - 20 mJy while 
the 87GB limit is ~ 25 mJy. Target sources were selected 
with a ^ —0.5 between 325 MHz and 4.85 GHz. There were 
670 CLASS-1 sources selected using WENSS (designated as 
CLASS-lb). The remaining CLASS-1 sources were selected 
on the basis of flux densities from the Texas 365 MHz cata¬ 
logue. Because the Texas survey has a flux density limit of 
250 mJy at 365 MHz, we relaxed the limits to 50 mJy or 


greater at 4.85 GHz and a ^ —0.6 between 325 MHz and 
4.85 GHz. All sources meeting the criteria were observed for 
& > 45° 20', and down to S ^ 35° south of the galactic plane. 
In all, 2550 sources were observed with this selection (des¬ 
ignated CLASS-la). Note that the Texas survey only covers 
the region S < 71?6. North of this declination, there was 
effectively no spectral selection in CLASS-la. 

In CLASS-2, we were able to select more sources against 
WENSS as that survey neared completion. In all, 4435 more 
targets were observed. The CLASS-2a session observed 2337 
targets that were nearly all selected using WENSS. The 
CLASS-2b observations mostly targeted sources below the 
WENSS southern limit (-1-30°) and were selected using the 
Texas survey and the White & Becker 1.4 GHz catalogue 
(which has a flux density limit of 100 mJy at 1.4 GHz), ex¬ 
cluding targets with a < —0.6 between 365 MHz or 1.4 GHz 
and 4.85 GHz. 

2.2 CLASS-3, CLASS-4, and the CLASS 
Complete Sample 

Our starting point for source selection in CLASS-1 and 
CLASS-2 had been the 87GB catalogue. However, the im¬ 
proved Green Bank survey GB6 became available in 1996. 
This catalogue contains 75162 sources in the declination 
range 0° ^ <5 ^ 75° with a 4.85 GHz flux density of 18 mJy 
or more and angular diameter 10^5 or smaller. Since the 
source positions and flux densities are better determined 
than in 87GB, GB6 was adopted as the primary finding 
list for CLASS-3, with the intention of bringing the entire 
CLASS survey into line with this new parent catalogue. In 
addition, the NVSS became available, giving coverage of 
the entire GB6 region at 1.4 GHz. Use of NVSS created 
a sample with a uniform rigorous spectral selection over vir¬ 
tually the whole northern sky. This GB6-targeted, NVSS 
spectrally-selected survey will be designated as the “CLASS- 
NVSS Complete Sample” or “complete sample”, and natu¬ 
rally updates the spectral selection used in JVAS as well as 
early CLASS. We used the NVSS version released as cata- 
LOG39.FIT for sample selection. This release purported to be 
complete in the GB6 region (catalog40.fit is now available 
filling in holes in other regions of the sky). 

The complete sample consists of 11685 targets with a 
GB6 flux density of 30 mJy or more, an NVSS source within 
70" of the GB6 position, and a spectral index a ^ —0.5 ver¬ 
sus NVSS. The spectral index is defined by adding the flux 
densities of all NVSS sources within 70" of the GB6 position 
and ignoring all other NVSS sources. The primary goals of 
CLASS-3 and CLASS-4 were to complete the sample. A list 
of 4982 sources for observation in CLASS-3 was drawn up 
from those sources not within 2' of a CLASS-1 or CLASS- 
2 pointing center, JVAS target, or identified with a known 
radio source (eg. from the VLA calibrator lists). 

In addition, the CLASS-3 list was supplemented with a 
list of GB6 sources (30 mJy or brighter), that have a WENSS 
source within 70" of the GB6 position, and spectral index 
a ^ —0.5 versus WENSS. This enables the construction of a 
“CLASS-WENSS Complete Sample” and will catch targets 
that have GHz-peaked spectra where the 1.4 GHz flux den¬ 
sity is above the spectral limit for inclusion in the CLASS- 
NVSS sample. 

In total, 5053 new sources were observed in CLASS3 
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Session 

Dates 

Yobs 

Ydet 

% detect 

Ymult 

% mult 

BW MHz 

Selection 

JVAS 

Feb 1990 - Dec 1992 

2720 

2613 

96.1% 

456 

17.5% 

25/50 

GB 

CLASS-la 

Feb 1994 - May 1994 

2550 

2199 

86.2% 

644 

29.3% 

25 

Texas 

CLASS-lb 

” 

670 

624 

93.1% 

65 

10.4% 

50 

WENSS 

CLASS-2a 

Aug 1995 - Sep 1995 

2337 

2150 

92.0% 

176 

8.2% 

50 

WENSS 

CLASS-2b 


2098 

1813 

86.4% 

230 

12.7% 

25 

Texas/WENSS 

CLASS-3a 

Feb 1998 - May 1998 

2338 

2043 

87.4% 

230 

11.3% 

50 

NVSS 

CLASS-3b 

” 

2058 

1823 

88.6% 

249 

13.7% 

50 

NVSS 

CLASS-3C 


657 

575 

87.5% 

125 

21.7% 

50 

NVSS 

CLASS-4 

Aug 1999 

1075 

753 

70.0% 

371 

49.3% 

50 

NVSS 

Total 


16503 

14593 

88.4% 

2546 

17.4% 



Complete Sample 


11685 

10906 

93.3% 

1558 

14.3% 




Table 2. JVAS/CLASS Statistics . These numbers reflect pointings that differ by more than 120". Duplicate pointings have been 
excluded. Note that a multiple source will be recorded in cases where bad data or mapping errors have resulted in a spurious secondary. 
The “% mult” column refers to the fraction of detected sources with multiple components of any flux density including these spurious 
multiples — if restricted to 2 mjy or brighter the total is 15.1% of sources detected for all sources and 11.6% in the complete sample. 


JVAS/CLASS Raw Component Separations 



Log Separation (arcsec) 


Figure 1. Histogram of maximum component separations for 
the 2219 CLASS sources that have at least 2 components with 
individual flux densities of 2 mJy or more (solid). Also shown 
are cuts at 10 mjy on component flux densities for 938 sources 
(dotted), and with component sizes of 120 mas or smaller for 333 
sources (dashed). 


from Feb - May 1998. The remaining sources that make up 
the complete sample were observed in CLASS-4, with a total 
of 1075 targets observed in August 1999. Note that most of 
the CLASS-4 targets were near the lower flux density limit 
of the sample, entered into the sample on the changeover to 
the NVSS, or were not detected in NVSS, and thus there 
are a higher percentage of non-detections and extended or 
multiple component sources (see § ^ and Table in this 
part of the survey. The final combined CLASS and JVAS 
archive contains 16503 distinct pointings (more than 120” 
apart), of which 11685 form the complete sample. The final 
complete sample covers an area of 4.96sr. 


3 CLASS VLA OBSERVATIONS AND 

ANALYSIS 

Two independent IFs of 25 MHz bandwidth (or 50 MHz 
with WENSS positions) were centered at 8.415 GHz and 
8.465 GHz (average 8.44 GHz) for GLASS-1 and CLASS-2. 
In CLASS-3, which was observed in its entirety using NVSS 
positions and 50 MHz bandwidth, the IF band centers were 
moved to 8.4351 GHz and 8.4851 GHz (average 8.46 GHz) in 
line with NRAO recommendations to avoid radio frequency 
interference — we will use 8.46 GHz as the fiducial CLASS 
frequency. For all CLASS observations, an on-source dwell 
time of 30 seconds was used, with 3.3-second integrations. A 
compact source from the JVAS calibrator list was observed 
after every n-th target source for phase calibration, where n 
was chosen based on the weather conditions prevalent during 
the session (for CLASS-1: n = 14; CLASS-2: n = 8; CLASS- 
3/4: n = 10-12). The JVAS calibration sources have rms 
position errors in each coordinate of 12 mas, so the result¬ 
ing CLASS source positions should have accuracies of about 
20 mas (in stable observing conditions) and hence should 
themselves be useful as potential phase calibration sources. 
We were able to observe one target source per minute, in¬ 
cluding the overhead from observing calibration sources and 
slewing between sources. 

The gain calibration using the JVAS astrometric sources 
was carried out in the Astronomical Image Processing Soft¬ 
ware (AIPS) package, distributed by NRAOQ The mapping, 
deconvolution and self-calibration were performed using the 
DIFMAP package (Shepherd, Pearson & Taylor 1994; Shep¬ 
herd 1997). The imaging of the sources was completed in 
an automated fashion using a DIFMAP script. The average 
rate for processing with the script was one source every two 
minutes. This procedure yielded images of good quality for 
most of the target fields. 

The relevant details of the calibration and mapping pro¬ 
cedures adopted for GLASS are discussed in the remainder 
of this section. 


^ http://www.cv.nrao.edu/aips 






















6 S. T. Myers et al. 


Statistical Sample - Flux Density 



Log Flux (mJy) 


Figure 2. Histogram of CLASS 8.46 GHz flux densities for detec¬ 
tions in the complete sample with components merged (solid) and 
GB6 4.85 GHz flux densities for the complete sample (dashed). 
The GB6 flux density limit of 30 mJy is easily seen in the dashed 
histogram. 


3.1 Calibration and positional accuracy 

Editing and calibration of the data was done using AIPS 
following the standard procedure. Initial processing was car¬ 
ried out at the time of the observations in order to generate 
prompt candidate lists for immediate follow-up where pos¬ 
sible. After the entire CLASS survey was completed, final 
consistent editing/calibration was done for the full dataset 
(including JVAS). The data were examined and edited us¬ 
ing AIPS TVFLG to remove obviously discrepant telescopes 
or baselines. The flux density scale was set (AIPS SETJY) 
using the standard calibrators 3C48 and 3C286, for which 
the flux density scale at 8.41-8.48 GHz in the observations 
was normalized to the Baars et al. (1977) values. Gare was 
taken to ensure that the flux densities on short baselines of 
these two sources approximated well to the established val¬ 
ues after the final correction was applied; corrections were 
applied by hand if necessary. 

The task AIPS GALIB was used to obtain telescope 
gain solutions on calibrators, and the resulting solutions 
were examined. Bad solutions were discarded by hand before 
interpolation by AIPS CLGAL into a final gain solution. Cal¬ 
ibration was applied to the CLASS sources and each CLASS 
source was split from the file into an output disk UV-FITS 
file. 

Fifteen sources were observed twice during the CLASS 
1-3 sessions. This provides a cross-check of the repeatabil¬ 
ity of our position determinations. These observations were 
calibrated independently, and did not necessarily use the 
same phase calibrators. For this set of sources, the rms off¬ 
set between the positions in the two epochs was found to be 
35.2 mas. This is slightly higher than, but consistent with, 
our expected uncertainties due to the JVAS position errors. 
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Figure 3. Histogram of spectral indices between CLASS and 
GB6 (solid) and between GB6 and NVSS (dashed). Sources from 
the complete sample detected in CLASS at 8.46 GHz are plotted 
(with any multiple components merged). The a ^ —0.5 selection 
between 1.4 and 4.85 GHz is clearly seen in the dashed histogram. 


3.2 The auto-mapping procedure 

The automap algorithm can be outlined as follows. An initial 
primary cycle is used to find peaks in a large image above 
some signal-to-noise ratio (SNR) cutoff. In a secondary cy¬ 
cle, small images are made around the location of each peak 
and a GLEAN window is created. In a tertiary cycle, decon¬ 
volution and self-calibration is carried out iteratively. The 
locations of the peaks (and thus the clean boxes) are stored 
for later use. 

A series of parameters define the procedure. The dimen¬ 
sions (cellsize and number of pixels) of the large primary im¬ 
age control the area around the pointing position in which 
components can be identified. In practice bandwidth smear¬ 
ing limits our effective field to about 1' for the 50 MHz ob¬ 
servations which form the majority of the sample, and thus 
a 2048 X 2048 pixel search image was constructed with 0'.'25 
cell size (giving a 256” x 256” searchable inner quarter) and 
super-uniform weighting. This field was searched for pixels 
with values above a given SNR cutoff above which sources 
are deemed to be reliably identified as real for further pro¬ 
cessing; SNR=8 was adopted for CLASS. 

Around each source, a secondary 2048 x 2048 image 
covering 140" x 140" (O'.'0684 pixels) was then made using 
uniform weighting, of which the inner quarter (70” x 70”) 
was cleaned. If a peak is found in this secondary image 
above a second SNR cutoff (the “field SNR cutoff”, for which 
SNR=10 was adopted), a CLEAN box with a size 75% of 
the fitted CLEAN beam size is placed around the image and 
25 iterations of the CLEAN algorithm are carried out with 
a gain of 0.05. The CLEAN component signal is subtracted 
from the data, and the data inverted back into the u-v plane, 
with the CLEAN procedure repeated until a residual cutoff 
level (5% of the peak or 6a) is reached. This means that 
we have a slightly higher sensitivity to secondary features 
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Table 3. First page of CLASS catalogue. Column 1 (bytes 0-16); GB6J name from the GB6 catalogue (Gregory et al. 1996); Column 
2 (18-25): Name in observing list; Column 3 (27-30); Reliability number of map (see text). NOSH indicates a non-detection; Column 
4 (32-42): Pointing position RA (J2000); Column 5 (46-55): Pointing position Dec (J2000); Column 6 (60-70): Observing epoch (90- 
92=JVAS, 94=CLASS1, 95=CLASS2, 98=CLASS3, 99=CLASS4); Column 7 (74-86): Detected component position(s) RA (J2000); 
Column 8 (89-101): Detected component position(s) Dec (J2000); Column 9 (103-111): 8.46-GHz flux density (mjy) of component(s); 
Golumn 10 (113-118): Fitted major axis (mas) of component(s); Column 11 (120-123): Fitted axis ratio of component(s); Column 12 
(125-128): Fitted position angle (degrees) of component. 


within this smaller held around the primary component than 
in the larger search held, due to the extreme taper applied 
to the visibility data input to the low-resolution search im¬ 
age (even accounting for the slightly lower SNR cutoff in 
the small field), and thus are somewhat more likely to find 
a faint source as a companion to a brighter source than on 
its own. 

For each field identified and CLEANed, a phase-only 
self-calibration with 10-second solution interval is performed 
and the CLEAN procedure is repeated. The CLEAN model 
is saved, and a model-fitting procedure (DIEMAP MOD- 
ELFIT) is run for 5 iterations on the u-v data. The start¬ 
ing point for model-fitting is the clean boxes found and 
stored during the iterative cleaning described above, with 
a Gaussian component assigned to each distinct CLEAN 
box (boxes must be further than 0"2 apart). The data prod¬ 
ucts of the auto-mapping procedure outlined above are the 
CLEAN model, the set of CLEAN windows, the list of MOD- 
ELFIT components, and the log-file generated by DIEMAP 
containing the commands executed plus output information 
from the noise tests (this log-file can be executed as a script 
in DIEMAP reproducing the original results). 

We have adopted the Gaussian MODELFIT results as 
the best representation of the DIFMAP analysis. Each com¬ 
ponent has a best-fit position, major axis size, elliptical 


Gaussian axis ratio and position angle. In addition, uncer¬ 
tainties for the fitted quantities are evaluated from the co- 
variance matrix. These quantities can then be filtered to 
select candidates with multiple components, and further fil¬ 
tered to find those with multiple compact components. Note 
that once distinct CLEAN boxes are laid down by the au¬ 
tomapper, a component is necessarily created that will be 
passed to the model-fitting. The model-fitting procedure 
may then move this component where it pleases, possibly 
inside the exclusion box of 0'.'2, or possibly outward to arbi¬ 
trarily large separations. There is no provision for deleting 
poorly constrained components during model-fitting, and 
these must be dealt with during the analysis phase. Eor ex¬ 
ample, a moderate flux density cutoff of 2 mJy per com¬ 
ponent, which corresponds to 5cr for good data taken with 
50 MHz bandwidth, is found to reject most of the extraneous 
components. 

The limitations of this algorithm include a restricted 
search field (only 256” square around the pointing center in 
the inner quarter of the large image), a relatively conserva¬ 
tive search SNR limit (8cr to avoid many spurious detections) 
in this large field, and dynamic range problems (often spuri¬ 
ous multiple components are found on sidelobes from bright 
sources). However, the automatic mapping has proved to be 
a reliable method of finding the best candidates quickly. 
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In some cases the automatic mapping program finds 
spurious structure or sources. Without manual intervention 
it is difficult to find the cause of many of these anomalies. 
However, the common symptom of such problems is that the 
final peak flux density after cleaning and self-calibration is 
much greater than the initial peak flux density. The ratio of 
these two numbers, the “reliability number” (RN), should 
theoretically be around 1. In practice we have found that 
maps with RN>1.5 usually indicate a data problem, and 
RN=2 almost certainly indicates a spurious detection. 

Because the goal of CLASS is to identify lens systems 
with multiple compact components, it is crucial that we eval¬ 
uate the eflectiveness of our automapper in finding such 
objects. To this end, we constructed fake visibility data 
sets by adding pairs of point components (using the AIPS 
task UVSUB) to actual CLASS fields in which no sources 
were detected by the VLA. This was repeated over a grid 
of component separations (0'.'15-10"), total flux densities 
(30-100 mjy) and flux density ratios (1:1-20:1) to simulate 
two-image gravitational lens systems. The pairs were ran¬ 
domly oriented on the sky and placed up to 20 arcseconds 
from the field center to mimic pointing inaccuracies. The 
data were then passed through the automapper, and the re¬ 
sulting MODELFIT output was analyzed. We find that for 
separations > 300 mas, the automapper consistently picks 
up two compact components when the fainter component 
has flux density > 2 mJy. For a source with 30 mJy of to¬ 
tal flux density, this corresponds to a flux density ratio of 
approximately 15:1. As expected, the identification of multi- 
component sources with larger total flux density is complete 
to larger flux density ratios. The automapper becomes in¬ 
creasingly inefficient at identifying multiple components as 
their separation decreases below 300 mas. At 150 mas sepa¬ 
ration, almost all of the data sets are fit to a single compo¬ 
nent, whose major axis increases with decreasing flux density 
ratio. 

Finally, there were found to be minor irregularities in 
the automapping for multiple component sources, with ex¬ 
traneous components being generated due to the high side- 
lobes of the VLA snapshot beam, particularly sidelobes from 
a source outside the small (70") field. This effect can be seen 
in the histogram of maximum component separations in the 
full sample is shown in Figure The secondary peak at 
around 20 arcseconds is almost certainly due to sidelobes, 
as the dirty beam for a VLA snapshot has its main side¬ 
lobes, which have an amplitude of 40% of the peak, in a 
hexagonal pattern spaced at 20" around the center. This 
was not deemed to be a serious problem for the lens sur¬ 
vey, as no true multiple sources would be lost, and simple 
filtering on the flux density and/or size of the modelfit com¬ 
ponents will remove most of these cases, as shown by the 
dotted and dashed curves in Figure |^. However, the reader 
must be warned that the automatically generated catalogue 
still contains these spurious sources. 


4 CLASS RESULTS AND STATISTICS 

In Table ^ we show the first page of the JVAS/CLASS results 
data base. The entire sample is available as an ASCII te xt file 


from [http : //www. jb.man. ac .uk/research/gravlens , This 
table contains the raw output from the automapping, and 


does not include manual remapping of questionable cases. 
Note that the last six columns contain the modelfitting in¬ 
formation for position, flux density, size and orientation of 
the Gaussian components. The axis ratio is often found to 
be zero, which indicates a component that it thinks is unre¬ 
solved in one dimension. If the major axis is small (150 mas 
or less), then this indicates an unresolved component, while 
if the major axis is significantly larger than this, then a zero 
or very small axis ratio usually indicates a poorly fitted com¬ 
ponent. These pathological cases are most often seen in cases 
where a large number of spurious and usually weak com¬ 
ponents are fit to uncleaned sidelobes on poorly calibrated 
data. For these objects, manual imaging, self-calibration and 
more careful modelfitting are being carried out. The reader 
should be careful in interpreting the auto-mapping output 
for sources with more than two components for this reason. 

Figure ^ shows a histogram of the 4.85 GHz flux den¬ 
sities of the complete sample taken from the GB6 survey. 
The 30 mJy lower flux density limit is clearly visible. Also 
shown in Figure ^ is the distribution of measured 8.46 GHz 
flux densities for detections in the complete sample. Sources 
generally have less detected flux density at 8.46 GHz than 
the 4.85 GHz survey flux density. This is for a combination 
of reasons but mainly spectral index and the fact that some 
of the sources are extended and partially resolved with the 
VLA in A-configuration. 

Histograms of the GB6 vs. NVSS spectral indices and 
CLASS vs. GB6 are plotted in Figure H for those sources 
in the complete sample that were detected in CLASS. 
The spectral index distribution derived using the observed 
8.46 GHz band merged flux densities and the 4.85 GHz GB6 
catalogue values is rather broad, centered around a « —0.6, 
with a long tail steepward. This behavior might reflect 
turnover in the spectra, source variability, resolution effects, 
a Malmquist-like bias, or likely a combination of most of the 
above. In particular, near the cutoff GB6 flux density the 
distribution will likely be dominated by sources scattered 
upward in flux density into the sample due to variability or 
measurement error over and above those lost downward out 
of the sample — this is the Malmquist bias effect. In addi¬ 
tion, steeper-spectrum sources tend to be extended and thus 
are more likely to be resolved in the high-resolution CLASS 
observations and to have depressed 8.46 GHz flux densities, 
further steepening their apparent integrated spectrum. 

In Figure ^ the measured 8.46 GHz flux densities are 
plotted versus the GB6 flux density for the complete sam¬ 
ple. The long tail toward steep spectral indices a < —0.5 and 
non-detections in GLASS are clearly shown. A similar plot 
of CLASS/GB6 (8.46 GHz to 4.85 GHz) versus GB6/NVSS 
(4.85 GHz to 1.4 GHz) spectral indices is shown in Fig¬ 
ure The spread in CLASS/GB6 spectral index is large. 
Furthermore, those sources with highly inverted (a > 1.5) 
spectra between GB6 and NVSS tend to have moderate 
GLASS/GB6 indices. The GLASS non-detections are again 
shown as a bar at the bottom, and these are predominantly 
sources with GB6/NVSS spectral index near the cutoff, al¬ 
though there is a significant tail of sources with inverted 
low-frequency spectra. The GLASS/GB6 ^ectral index is 
plotted versus GB6 flux density in Figure [^ The trend to¬ 
ward larger scatter in spectral index with lower GB6 flux 
density is clearly seen. The upper envelope with inverted 
spectra concentrated to lower flux densities at the lower fre- 
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Figure 4. Observed JVAS/CLASS 8.46 GHz flux density versus GB6 4.85 GHz flux density, for the complete sample. For sources found 
with multiple components, the merged 8.46 GHz flux densities were used. Sources above the dashed line have apparent spectral indices 
between CLASS and GB6 flatter than a > —0.5. The bar of sources at the bottom represents GB6 targets where no CLASS source was 
detected (the assigned 8.46 GHz flux density is arbitrary and randomly spread to show the distribution with 4.85 GHz flux density). 


quency band is expected, however the excess of steep spec¬ 
trum sources at the lower flux densities is worrisome, and 
likely related to the problem of non-detections. Finally, the 
8.46 GHz to 1.4 GHz NVSS spectral index is plotted versus 
the NVSS flux density in Figure The spectral index range 
is tightened up considerably, and the envelopes due to selec¬ 
tion are clearly seen with the non-detections clustering at the 
limit for the spectral and flux density cutoffs. The “cleaner” 
nature of this relation is likely due to being well above the 
NVSS survey flux density limit (reducing bias), the fact that 
the GLASS and NVSS surveys are nearly contemporaneous, 
and the intrinsic reduction in variability effects at the lower 
NVSS frequency. 

There is a known increase in the rms uncertainty in 
GB6 flux densities at lower declination, with the 5a limit So 
increasing from less than 20 mJy at declination S = -1-20° to 
over 30 mJy at the equator (Gregory et al. 1996). To investi¬ 
gate this effect, we plot the histograms of the GLASS/GB6 
spectral index splitting the sample at declination -1-20° (Fig¬ 


ure H). There is a slight widening of the distribution and 
shift of the centroid toward steeper spectral index in the 
lower declination sub-sample consistent with the effect. This 
points toward problems with GB6 flux densities near the se¬ 
lection cutoff, especially as we see a significantly increased 
non-detection rate at lower declinations also (see below). 

4.1 Detection rate 

The statistics of detections and multiple components found 
in the JVAS/GLASS sessions are given in Table Of the 
16503 targets observed, 14593 were detected in the auto¬ 
mapping stage (88.4%). 

In Figure 0 we display the percentage of observed 
sources which were not detected in the complete sample, 
as a function of 5-GHz flux density and spectral index. It 
can be clearly seen that the detection rate drops both with 
decreasing flux density and steeper spectra, as expected. In 
the survey as a whole, the non-detection rate is much higher 
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Figure 5. Spectral indices between CLASS and GB6 (8.46 GHz and 4.85 GHz) plotted versus GB6 to NVSS (4.85 GHz and 1.4 GHz) 
index for the complete sample. As before, multiple components are merged. The bar of open circles at the bottom represents sources in 
the sample for which no 8.46 GHz detection was made, plotted versus the selection spectral index (y-values randomly spread to show 
the distribution). 


(11.6%) than in the complete sample (6.7%). The extra non¬ 
detections arise from a variety of reasons which have already 
been described. The main extra contribution is likely to arise 
from parts of the survey where the spectral selection was 
relaxed or absent (particularly parts of CLASS-1). This is 
likely to have caused significantly more contamination by 
steep-spectrum and resolved sources. In addition, a few of 
the inverted-spectrum non-detections found at lower galac¬ 
tic latitudes are identified as planetary nebulae. 

The 6.7% non-detection rate in the complete sample is 
relatively low: nonetheless a 30-mJy source at 5 GHz should 
be detectable at > 20a at 8.46 GHz, even with a moderately 
steep radio spectrum. If we split the complete sample at 
declination -1-20° as we did earlier (see Figure we find that 
for S > 20° the non-detection rate is only 3.8%, and rises 
to 10.4% for 5 < 20°. Again, this points toward problems 
related to the GB6 flux density scale near the cutoff, as 
the GB6 rms flux density uncertainty rises sharply at low 
declinations. Note that the non-detection rate for GB6 flux 


densities above 60 mjy drops to 2% overall (46% of the total 
sample is above this cutoff). 

Of the non-detections, 100 are within the existing re¬ 
gions of the FIRST survey (Becker, White & Helfand 1995) 
and can be investigated in detail. Just over half of the non¬ 
detections (53/100) are due to classical double radio sources, 
in which one of the radio lobes is not picked up within the 
70” search radius in the GB6/NVSS correlation (if the sec¬ 
ond lobe is included, then the source would be dropped from 
the sample as the combined spectral index will pass the cut¬ 
off). Of the remaining 47 sources, 21 are revealed as classi¬ 
cal double sources by FIRST. These are either single NVSS 
identifications, in the case of doubles with separation < 45”, 
or more rarely larger sources; in either case they appear in 
the GLASS survey because the overall spectral index falls 
just under the a\'T' cut. 26 sources remain. Of these, eight 
are obviously large lobes or extended structure in which in¬ 
dividual hotspots have registered in NVSS, resulting in in¬ 
clusion in the CLASS sample. For eighteen sources there is 
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Figure 6. Spectral indices between CLASS and GB6 (8.46 GHz and 4.85 GHz) plotted versus GB6 4.85 GHz flux density for the 
complete sample. As before, the bar of of open circles at the bottom represents sources in the sample for which no 8.46 GHz detection 
was made. 


no obvious explanation for why they were not detected at 
8.46 GHz; they generally have spectra approaching the spec¬ 
tral index cut, but are unresolved or only slightly resolved 
by the FIRST survey, and 6 of these 18 sources have a sin¬ 
gle NVSS identification and are unresolved by FIRST. How¬ 
ever, these unexplained non-detections are most likely to be 
either very small, relatively flat-spectrum double sources, 
highly variable flat-spectrum sources, or sources for which 
the observations failed (although there is no evidence for the 
latter). In general, therefore, we can state that the contami¬ 
nation of the survey by intrinsically steep-spectrum sources 
is about 4-5%, and about 1.5% of the survey sources are not 
detected for reasons that are not clear. 


4.2 Statistics of the northern CLASS-1 region 


As mentioned in §2.1, the Texas 325 MHz survey only cov¬ 
ered the region S < 71? 6 and thus north of this declina¬ 
tion there was effectively no spectral selection in CLASS-la. 
There were 489 87GB-selected sources that were observed in 


this manner, out of a total of 533 sources in 87GB in this 
region and flux density of 50 mJy or higher, with the re¬ 
maining 44 sources having been previously detected in other 
surveys. These 533 sources form a complete sample free of 
spectral selection that can be used to check the efficiency of 
our selection criteria in filtering out unresolved or complex 
sources. 


Of the 489 CLASS targets in the north region, 350 were 
detected in the automatic mapping. Added to the 44 pre¬ 
vious detections in the region, the total detection rate is 
394/533 (73.9%). This detection rate is significantly lower 
than the average of the entire CLASS survey (89%). In ad¬ 
dition, 155 of the 350 CLASS detected sources were found 
to have multiple components (44.3%), a significantly higher 
rate than the overall CLASS average (14.3% of detected 
sources in the complete sample). These statistics are consis¬ 
tent with the expectations for the spectral selection, which is 
designed to reduce the number of extended and intrinsically 
complex sources in the sample. 
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Figure 7. Spectral indices between CLASS and NVSS (8.46 GHz and 1.4 GHz) plotted versus NVSS 1.4 GHz flux density for the 
complete sample. As before, the bar of of open circles at the bottom represents sources in the sample for which no 8.46 GHz detection 
was made. 


5 DISCUSSION AND SUMMARY 

In total, 15.4% of the entire survey and 13.3% of the com¬ 
plete sample (17.4% and 14.3% of detected sources respec¬ 
tively) have multiple components identified in the automap¬ 
ping process. Some of these multiple sources, however, will 
be due to sidelobes and other artifacts of the automapping 
process. For further analysis we require an estimate of which 
sources are genuine multiple objects. For lens candidate se¬ 
lection, in particular, further constraints are applied in total 
8.46-GHz flux density, in primary-to-seconda^ flux density 
ratio and in component separation. In Figure fy we show the 
histogram of maximum component separations in the full 
sample (the complete sample shows essentially the same dis¬ 
tribution). The survey is clearly sensitive to multiple compo¬ 
nents with separations from 160 mas to 300 arcseconds, and 
thus our primary goal of identihcation of gravitational lens 
candidates is achievable. Note that in Figure double-lobe 
sources, knots within lobes, core-jet sources, superposition 
of unrelated sources, artifacts of the automapping process 


and real multiple-image lens systems are all mixed together, 
and lenses are a tiny fraction of the total number of multi¬ 
ples! When restricted to components with flux densities of 
2 mJy and above, 15.1% of all sources detected and 11.6% 
of detected sources in the complete sample were multiple — 
these numbers better represent the fraction of true multiples 
in the survey. Details of candidate selection will be discussed 
in Paper 2 of this series. 

Optical identihcations have been sought from scans of 
the Palomar Sky Survey prints carried out by Kibblewhite 
et al. (1984). As expected from earlier work, the identih- 
cation fraction of fiat-spectrum sources to R~20 decreases 
with decreasing flux density, reaching 50% at ~ 50 mJy 
and decreasing further thereafter. Redshift surveys are in 
progress to identify sources at these flux density levels 
(Falco, Kochanek & Munoz 1998; Marlow et al. 2000). Fig¬ 
ure ^ shows the optical identihcation rate to this level; it 
can be seen that the decrease with 5-GHz flux density is 
monotonic, and does not depend strongly on radio spectral 
index. 









CLASS: source selection and observations 13 


Statistical Sample - Spectral Index 



Spectral Index a(4.85 —8.46) GHz 


Figure 8. Histogram of spectral indices observed between CLASS 
and GB6 as a function of declination: Dec > 20° (solid) and Dec 
< 20° (dashed). The vertical axis is the fraction of the total falling 
in that bin. Only sources from the complete sample with CLASS 
detections are plotted (with any multiple components merged). 

In summary, we have presented radio observations of 
a complete sample of flat-spectrum radio sources with the 
VLA. A total of 16503 distinct target sources have been 
observed, of which 11685 form a complete sample with the 
properties Ssghz > 30mJy, |fe| > 10°, > —0.5. Con¬ 

tamination by steep-spectrum sources is around 4% in the 
complete part of the survey. Further papers in the series will 
discuss the selection of candidate gravitational lenses and 
astrophysical inferences from the statistics of this sample. 
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